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THE  DETERMINATION  OF  LATENT  ROOTS  AND  INVARIANT  MANIFOLDS 
OF  MATRICES  BY  MEANS  OF  ITERATIONS1 
K.  A.  Semendiaev  (Moscow) 


In  the  practical  solution  of  many  problems  of  vibration  theory 
the  greatest  difficulty  is  occasioned  by  the  necessity  of  solving 
the  so-called  characteristic  equation.  Among  many  methods  long  since 
suggested  for  the  solution  of  this  and  other  associated  problems, 
one  of  the  most  effective  and  widely  used  is  the  method  of  iteration. 

Given  that  this  method  is  to  be  utilized,  the  presence  in  the 
matrix  in  hand  of  characteristic  roots  moduli  of  which  are  equal, 
as  well  as  the  character  of  the  invariant  manifolds  of  the  matrix, 
substantially  influence  the  course  of  the  process  of  successive 
approximations.  In  works  devoted  to  the  method  of  iteration,  a 
full  analysis  has  not,  up  to  the  present,  been  adduced  of  all 
cases  possible  under  the  circumstances. 

^Prikladnaia  matematika  i~mekhanika  (Institute  of  Mechanics,  USSR 
Acadeny  of  Sciences,  Applied  Mathematics  and  Mechanics),  vol.  7 
(19U3),  pp.  293-322  (original  in  Russian,  with  English  summary). 
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R.  Kises  and  H.  Pollaczek-Geiringer ^ ^ ^ limit  thx. selves  to 

symmetric  matrices  having  no  multiple  latent  roo'os.  D mean  and 

[2  31 

Collar  * J admit  the  presence  of  a single  pair  of  conjugate  com- 
plex latent  roots  of  maximum  modulus^  or  of  two  latent  roots  of 
coincident  moduli.  The  latter  case  is  considered  by  passing  to 
the  limit,  in  connection  with  which  the  possibility  of  a coinci- 
dence of  latent  roots  without  an  elevation  of  the  degree  of  the 
corresponding  elementary  divisors  is  not  taken  into  consideration. 

The  method  of  iteration  is  much  more  thoroughly  investigated 
in  a work  of  Aitken^,  which  contains  in  part  results  identical 
’.rith  those  of  f 5 and  7 of  the  present  work  • 

In  the  first  part  of  the  present  work  ( $ 1-9)  a complete 
analysis  is  carried  through  of  all  possible  cases  of  the  course  of 
iterations,  and  instr actions  are  given  for  practical  methods  of 
determ? r>iny  *-he  dominant  latent  roots,  as  wall  as  the  correspond- 
ing invariant,  ...  i. folds  of  the  matrix,  for  each  distinct  case. 
Essentially  the  investigation  is  carried  out  by  the  geometrical 
method*  A familiarity  with  the  fundamentals  of  the  theory  of 
elementary  divisors  and  of  n-dinensional  Euclidean  geometry  is 
required  of  the  reader.  (See,  e,gf ) 

The  second  part  of  the  work  gives,  for  the  most  general  case, 
a method  for  a successive  determination  of  the  latent  roots,  in 

(Translator’s  note)  "The  latent  root/s  of  greatest  modulus" 

(Russian)  is/are  spoken  of  as  "the  dominant  latent  root/s" 

(English). 

2 

Aitken’s  work  became  known  to  the  author  only  after  completion 
of  the  present  study. 
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order  of  diminishing  mod  ill,  and  of  the  corresponding  invariant 

manifolds.  The  procedure  is  a generalization  of  the  method  pro- 

[2  31 

posed  for  the  simplest  case  by  Duncan  and  Collar  # J and  also, 
independently  of  them,  by  P.  F.  Papkovich^* 

§ 1.  BASIC  CONCEPTS 
let  there  be  given  a matrix 

D * {dik}  * i*k  ■ 1»2»  — * n 

with  real  elements.  In  n-dimensional  Euclidean  space  with  a given 
orthogonal  basis  this  matrix  defines  a linear  transformation  that 
relates  to  any  vector  a(ap  , aR)  the  vector  b • Da,  the  co- 
ordinates of  which,  b^,  •••  , b^  are  defined  by  the  formulae 


b. 

i 


(i.i) 


The  vector  f ( e<,r-  o<  r,  •••  . * r)  is  called  the  latent 

-r  l * c w w n 

vector  corresponding  to  the  latent  root  Ar,  if 


Df, 


Arfr 


(1.2) 


Since  the  vector  is  by  this  equation  defined  only  up  to  a 

constant  multiplier,  we  shall,  unless  it  is  otherwise  specified, 

henceforth  consider  the  vector  f normalized  so  that 

~r 


^It  should  be  noted  that  the  conclusions  of  Papkovich 1 3 work  are 
suitable  only  to  a symmetric  matrix  whose  latent  roots  are  without 
coincident  moduli,  rather  than  to-  an  arbitrary  matrix  vrith  roal 
latent  roots  as  the  author  states. 
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The  latent 

roots  of 

a matrix  are 

characteristic 

equation 

dn  - A 

d12 

• • • 

dln 

d21 

d22  - X 

• •• 

d2n 

• • • 

• • • 

• • • 

• • • 

dnl 

dn2 

• • • 

d - X 
nn 

|D  - XE|  . 0 (1.3) 


where  E designates  the  diagonal  init  matrix.  To  simplify  the  ex- 
position we  shall  consider  the  matrix  D to  be  non^singular,  i.e., 
not  having  null  latent  roots  (|d|  ^ 0)«  Equation  (1*3)  has  n roots, 
Xl#  X2#  •••  # Xn#  which  we  may  enumerate  in  order  of  descending 

moduli,  i«e», 

| X1|  ^ | X2I  A I X3I  • •••  * I *nl 

Among  these  there  may  of  course  be  equal  roots* 

If  a latent  root  Ar  be  known,  then  the  determination  of  the 
corresponding  latent  vector  reduces  to  the  solution  of  the  system 
of  linear  homogeneous  equations 

(dik  " Xr  *i|c>*kP  “ 0 (i  ■ 1»  2#  •••  » n)  (1.U) 


^ik  ^onecker'3  delta. 

We  introduce  into  the  subsequent  discussion,  in  addition  to 
the  latent  vector s9  the  invariant  manifolds  of  a matrix*  A linear^ 

^Since  only  linear  manifolds  are  considered  in  this  work,  the  word 
linear  i/ill  in  future  be  dropped* 


■ 
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manifold  of  k (k  < n)  dimensions,  E^  lyirg  in  is  called  invari- 
ant  if  a vector  transformed  by  means  of  the  matrix  D from  any 
vector  belonging  to  also  belongs  to  E^* 

The  simplest  exangle  of  an  invariant  manifold  is  a straight 
line  whose  direction  coincides  with  that  of  a latent  vector* 

The  structure  of  the  invariant  manifolds  of  a matrix  in  real 
space  2^  essentially  depends  on  the  elementary  divisors  of  the 
matrix  D,  as  also  upon  the  presence  of  complex  roots. 


§ 2.  THE  INVARIANT  HAITIFOIDS  HI  THE  CASE  OF  REAL  LATE1IT  ROOTS 
Let  us  first  suppose  that  all  the  latent  Hoots  are  real. 

Here,  as  is  known,  the  following  cases  are  possible, 

1%  All  the  latent  roots  are  distinct;  to  each  latent  root 
corresponds  a determinate  latent  vector.  Any  linear  manifold  E^ 
constructed  on  k latent  vectors  is  invariant;  conversely,  any 
k-dinensional  invariant  manifold  contains  k independent  latent 
vectors  which  nay  be  taken  as  a basis  for  this  manifold, 

2,  There  ire  jiiultiple  latent  roots,  but  the  elementary 
divisors^"  of  the  matrix  D — XE  ai*e  relatively  prime.  To  each 
latent  root  A.^  of  multiplicity  r^  there  corresponds  an  r^-riinen- 
sional  invariant  manifold  Ep.  containing  a unique  invariant  direc- 
tion. For  a basis  of  this  Ep^  there  may  be  taken  the  system  of 
vectors  •••  , having  the  fo  Hot  ring  properties: 

Dfi1-  Vi1,  Dfj1-  1»  X- 2,3,  ...  ,r±  (2.1) 

Tn  this  work  ths  treatment  of  elomontary  divisors  follows 
Weierstrass, 


. 

■ 


6 


We  may  call  these  vectors  the  fundamental  basis  ( ocho$hom 

of  the  invariant  manifold  Er  that  corresponds  to  the  latent  root 


A 


i* 


For  simple  latent  roots  the  fundamental  basis  for  each  latent 
root  will  consist  of  a unique  vector. 

All  the  fundamental  bases  corresponding  to  the  distinct  latent 
roots  contain  in  the  aggregate  n vectors,  which  are  linearly  inde- 
pendent and  may  serve  as  a basis  of  E^. 

From  equality  (2.1)  it  follows  that  a linear  manifold  E^  con- 
structed from  the  vectors  •••  , (%  * ri)  will  be  an 

invariant  manifold. 

In  this  way  we  may  have  here  a system  of  invariant  manifolds 
lying  one  within  the  other,  of  two,  three  and  so  forth  on  up  to  r^ 
dimensions,  containing  a unique  invariant  direction  f^1.  Invari- 
ant manifolds  of  this  type  we  shall  call  axial  (oceffoe)  and  f^ 
the  axis  (ocb)  of  this  manifold. 

Apart  from  axial  invariant  manifolds,  in  this  case  it  is  pos- 
sible to  construct  invariant  manifolds  taking  as  a basis  an  arbi- 
trary number  of  each  of  the  first  vectors  of  the  fundamental  basis 
corresponding  to  the  distinct  latent  roots.  With  this  all  invari- 
ant manifolds  of  the  matrix  D have  been  exhausted. 

3*  The  elementary  divisors  of  the  matrix  D - AE  have  common 


divisors. 

For  example,  let  there  be  for  a certain  latent  root  the 
corresponding  elementary  divisors 

(A  - Aj/1,  (A  - Ah)*2,  •••  , (A  - Ah)6k,  ‘ «2  * — * *k 
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We  shall  call 


i-1 


The  n unber  of  elementary  divisors  with  exponent  j we  shall 
designate  by  itk  (1  * j * p)  and 

^ n.  « ky  , % " 1»  2,  •••  , p , k.  » k 


ITobe  that  of  course  not  all  m.  are  different  from  zero* 

0 

Then  to  the  latent  root  there  corresponds  a s^-dimen- 
sional  invariant  manifold  Es^  with  structure  as  follows:  there 
is  a k^-dimensional  manifold  in  which  each  direction  is 
invariant#  In  E^  lies  a k2**dinensional  manifold  Ej^  such  that 
through  each  direction  in  it  passes  a two-dimensional  invariant 
plane  not  belonging  to  Ej^# 

Iloreover,  to  each  mmber  of  domensions  % * p there  corresponds 
a manifold  of  invariant  directions  E^CIE^  such  that  each  direc- 
tion in  it  is  an  axis  of  determinate  axial  invariant  manifolds 
placed  one  within  the  other,  these  manifolds  being  of  2,  3 and  so 
forth  up  to  % dimensions  and  intersecting  E^^  only  along  the  axis. 

In  the  light  of  this,  obviously  Cl  if  % > j.  All  these 

A 0 

invariant  manifolds,  ending  with  the  p-dimensional,  completely 
fill  up  a linear  manifold  of  dimensions,  since  2^^  jm^  • s^# 
Any  vector  in  this  E^  belongs  to  a certain  axial  invariant  mani- 
fold of  dimension  not  greater  than  p.  As  a basis  of  this  manifold 
it  is  always  possible  to  select  a system  of  independent  vectors 


. 


V 
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satisfying  equations  analogo  is  to  equations  (2,1).  Such  a system 

we  shall  also  call  a fundamental  basis  of  the  given  axial  manifold. 

In  the  simplest  case,  when  * •••  ■ * 1,  we  shall 

have  a k-dinensional  manifold  E^  in  which  each  direction  will  be 

invariant.  Here  E^  ® E . 

K1  sh 

In  the  general  case  to  each  latent  root  will  correspond 

an  invariant  manifold  E (Z  s,  • n),  of  the  structure  indicated 
n h n 

above. 

An  arbitrary  invariant  manifold  of  the  matrix  D may  be  con- 
stricted upon  the  independent  axial  invariant  manifolds^,  the 
rcimber  of  dimensions  of  which  does  not  exceed  the  greatest  exponent 
of  the  power  of  the  elementary  divisors  of  the  matrix  D - X E be* 

longing  to  the  distinct  E_  . 

*h 

§3.  THE  INVARIANT  MANIFOIDS  IN  THE  CASE  OF  COMPIEX  LATENT  ROOTS 

Assume  now  that  there  are  complex  roots  among  the  latent  roots 
of  the  matrix  D,  Since  we  have  limited  ourselves  to  the  considera- 
tion of  matrices  with  real  elements,  the  complex  latent  roots  must 
necessarily  be  conjugate  pairs.  All  results  adduced  in  f 2 are 
transferable  to  this  case  in  complex  space,  but  for  real  space 
corrections  must  be  introduced. 

To  each  pair  of  conjugate  latent  roots  there  corresponds  a pair 
of  conjugate  complex  latent  vectors,  determining  a real  invariant 
two-dimensional  plane.  Thus  in  the  case  in  which  all  the  latent 
roots  are  distinct,  to  the  invariant  directions  (for  the  real  latent 
roots)  such  invariant  planes  are  added  in  the  role  of  fundamental 
elements  from  which  all  the  invariant  manifolds  are  construct able. 
^Several  independent  axial  manifolds  can  be  taken  from  each  E^. 
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In  the  second  case  (latent  roots  multiple,  elementary  divisors 
relatively  prime),  to  the  tx*o  conjugate  latent  roots  fy  exp  ( ±i0^) 
of  multiplicity  r^  there  correspond  invariant  manifolds  of  tx/o, 
four  and  so  forth,  up  to  2r^  dimensions,  located  one  xiithin  another. 

These  manifolds  may  be  constructed  by  means  of  a system  of  linearly 

k k k 

independent,  fundamental  vectors  (real),  f ^ , f 2 , •••  , fg£  , 

k 

satisfying  the  equations 

Df£-,k(cos  4^-sin  44' )>  DVlVk(cos  ^-l“Sin 

(3.1) 

Vk(sin  44+ cos  44),  m2%  Vk(si*  44^1 + «=os  44^)+4x-2 


t - 2,  3,  •••  , rk 


Invariant  manifolds  of  this  kind,  constructed  upon  two  conplex 

axial  manifolds,  ue  shall  call  imaginary-axial  («Hi«iO-oce»oe)  invari- 

k k k 

ant  manifolds,  and  the  system  of  vectors  f^  , f ^ 9 •••  > f^ 

k 

fundamental  basis  in  E2rk, 

The  structure  of  the  fundamental  invariant  manifolds  that  cor- 
respond to  the  complex  latent  roots  in  the  third  case  considered  in 
§ 2 xd.ll  change  accordingly.  We  shall  not  dwell  in  detail  upon 
this, 

h.  CERTAIN  PROPERTIES  OF  THE  SIMPIEST  INVARIANT  MANIFOLDS 

In  the  following  we  shall  consider  vectors  obtained  from  a 
certain  initial  vector  qQ  by  successive  transformations  effected 
by  the  matrix  D,  i.e,,  the  vectors 

3k  ’ Dk3o  » k - 3-,  2,  3,  ••• 


(U.D 


■ 
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Assume  first  that  the  vector  qQ  lies  in  a certain  imensional 
axial  invariant  manifold  E^  corresponding  to  the  latent  root  A^. 

Let  qQ  be  analyzed  according  to  the  vectors  f^1,  •••  , f^  forming 
the  fundamental  basis  of  E^: 

5o  * cl~l  + c2»2*  + **•  + (^»2) 

Here  it  may  be  assumed  that  c^  / 0,  i.e.,  that  qQ  cannot  be  placed 
in  an  axial  manifold  of  dimension  less  than  % in  number.  Then  on 
the  basis  of  equality  (2.1)  there  is  readily  obtained: 

Sp-  AiP(c1f1i+c2f2i  + ...+c^i)+  A.P“1(P)(o2f2i+c3f3i  + ...+c^1) 

+ Af2(P)(c3f3i  + •••  + c^f^)  + •••  + (U.3) 

From  this  equality  we  see  that  with  increasing  p, 

LV 


lim 


(U.U) 


and  accordingly  the  vector  q^  will  approach  f^1  in  direction . 

In  the  simplest  case,  when  qQ  is  collinear  with  the  latent 
vector  (X  * l)f  equation  (U.3)  takes  the  form 


Sp  • Ai  S' 


(U.S) 


It  is  known"*"  that  every  vector  q belonging  to  Ey  satisfies 


the  equation 


0>  - - 0 


(U.6) 


^ec,  c.g.,f6l,  p4  179, 
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i#e.,  that  the  vectors  qQ,  •••  , are  connected  by  a linear 
dependence 


whence  we  obtain 


or 


(r+1)  VSof'i-AW'S*1  “ (*"r) 1 3c3l‘ • • VlSr+lW • 

v ( 

r • 0,  1,  •••  , 


(4.9) 


Here,  as  in  what  follows,  [q^,  •••  q^]  designates  the  poly- 
vector constructed  upon  the  vectors  q^,  qg,  •••,q^.  [See  Note  1 
at  end  of  translation. ] 

We  now  adduce  analogous  formulae  for  Imaginary-axial  invariant 
manifolds.  Let  us  consider  first  the  simplest  case,  in  which  qQ 
lies  in  the  tw»«dimensional  imaginary-axial  manifold  that  cor- 
respo 


(4.10) 
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From  equations  (3.1)  we  obtain 

DpqQ  - AfkP[cos(*  +P0k)j‘1k  + sin(  + p0k)f2k]  (U.ll) 

From  (lull),  or  directly  from  (3#l)  and  (U*10)  it  follows  that 
(D2  - Zfk  cos  0kD  ♦ fk2)qo  . 0 (U.12) 

Whence  it  emerges  that 

[3ia2l  * <'k2[3o3i]*  [3o92]  ‘ 2^k  cos  4[3o3i]  (U‘13) 

Let  us  now  pass  to  the  general  case.  Let  qQ  belong  to  the 
imaginary- axial  invariant  manifold  "that  corresponds  to  the 
latent  roots  ^ exp  (±  i0k),  and  let  the  resolution  of  qQ,  in 
terms  of  the  vectors  forming  the  fundamental  basis  in  3,,^  have 
the  form  (we  drop  the  upper  indices  k): 

qQ*A^(cos  *]£i+sin  * + A2(cos  o^f^+sin  oc2f^)  + **»  + 

(U.1W 

+ A^(cos  o<  ^2^—1  + sin 
Then  on  the  basis  of  (3*1)  we  have 

Sp  -rkp  i A^.[cos(o<^  + p0k)f2^1  + sin(*  j + p0k)f2j]  + (U*15) 

+ (?vr  XAj+itcos(*‘j+rfp-lK)i2j*3+sin(*,j4i+{p"1K)^23-2]+',*'f 

♦ ^2^-l^A^k"^cos^X+  ♦ sin(*^*  {p-X+l}0k)f2] 
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Equation  (U.12)  may  here  be  replaced  by  the  equation 

(D2  - 2fk  cos  j^D  ♦ (»k2)^qo  - 0 (H.16) 

which  expresses  the  sought  relation  of  linear  dependence  between 
the  vectors  qQ,  qp  , q2^»  From  (U#l6)  it  follows  that 

C3o3l  •••  5i-l3r+l  S2J 0 * crf2X"r[3o3i  — 32Z-1]  (u*17) 


where 


'2m 


'2m+l 


so 


22i  cos21*1^ 


c„  - , | 22i+1cos2U1t 


(4.18) 


m*i  J \ 2i+l, 


the  summation  proceeding,  according  to  i,  from  0 to  min  (n,  /-m) 
in  the  first  sum,  and  to  min  (m,  X“m“l)  in  the  second. 

In  particular 


[3i32  •**  32X]  “ Pk^SoSi  32X-i] 


(4.19) 


From  formulae  (U*18)  it  is  readily  seen  that  c^  **  c2%^r* 

55#  THE  METHOD  OF  ITERATION  FOR  THE  DETERMINATION 
OF  THE  DOMINANT  LATENT  ROOT 

Let  us  now  designate  by  Xp  *2*  •••  9 Am  the  distinct 
latent  roots  of  the  given  matrix  D ■ {d^  J j further, 

| XjJ  * | X2|  4 ...  * | xj|  and  by  s^9  •••  , respectively 

their  multiplicity  ( Z s^  • n).  Let  the  rnmber  of  complex  latent 
roots  (distinct)  be  equal  to  2t. 


To  each  real  latent  root,  and  to  each  pair  of  conjugate  complex 
latent  roots,  there  corresponds  a proper  invariant  manifold,  as  has 
already  been  shown  in  1 2 and  3. 

Ho  two  of  these  manifolds  have  common  directions,  and  all  of 
them  together  determine  an  n-dimensional  space,  so  that  an  arbitrary 
vector  a may  be  resolved  into  components 

5 * *i  + + •**  + ~am  (5.1) 


such  that  each  component  belongs  to  a certain  axial  invariant 
manifold  Ep<(r^  * s^),  the  latter  corresponding  to  the  latent  root 
X If  a limitation  to  real  space  be  made,  for  the  complex  con- 
jugate latent  roots  X^  and  X^+^  « must  be  considered  as 

belonging  to  the  corresponding  imaginary-axial  invariant  manifold 
E2r^(r^  < s^),  whereas  • 0.  Depending  upon  the  selection  of 
the  vector  a,  separate  sunnands  of  sum  (5.1)  may  reduce  to  zero. 

Transforming  the  vector  a p times  by  means  of  the  matrix  D, 
we  shall  obtain 


Vpa  « tP^  + iP^  + ...  + Dpam 


(5.2) 


Each  summand  in  this  sum  can  be  expressed  in  terras  of  the  basic 
vectors  of  the  corresponding  invariant  manifolds  in  accordance 
with  formulae  (U.3)  and  (U*15).  From  these  formulae  it  follows 


directly  that  the  ratio 


ip^i 


will  tend  to  zero  if  | X^l  < | X^|, 


and  will  grow  — or,  in  the  simplest  case  (a^  coincides  with  the 
invariant  direction),  remain  constant  with  | XjJ  ■ | X^|.  [See 
Note  2 at  end  of  translation.] 


- 

- 
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Therefore  no  matter  what  the  assigned  degree  of  accuracy  nay 
be,  given  a sufficiently  large  p it  is  possible  to  consider  the 
vector  D^a  as  lying  in  the  invariant  manifold  formed  by  the  funda- 
mental axial  and  imaginary-axial  manifolds  that  correspond  to  the 
dominant  latent  roots,  if  only  the  vector  a contains  a counterpart 
in  this  manifold  (fundamental  theorem). 

If  in  sum  (5.1)  the  first  summand  different  from  zero  be  a^, 
then  A-^  should  be  replaced  by  A^  everywhere  above.  It  is  never- 
theless practically  always  possible  to  consider  a^  / 0.  Let  alone 
the  fact  that  the  probability  of  the  absence  of  a component  of  the 
arbitrarily  selected  initial  vector  in  some  determinate  invariant 
manifold  or  another  is  equal  to  zero,  it  should  be  pointed  out 
that  inasmuch  as  the  computations  are  being  conducted  to  a limited 
number  of  decimal  places,  the  component  a^  will  be  present  if  only 
because  of  the  rounding  error  — and  this  will  be  so  whether  or  not 
the  attempt  be  made  to  select  vector  a such  that  it  is  known  to 
contain  no  component^  a^.  This  circumstance,  as  will  be  shown  in 
what  follows,  considerably  complicates  the  determination  of  the 
latent  roots  whose  moduli  are  less  than  that  of  A 

The  method  for  the  approximate  determination  of  the  dominant 
latent  roots  by  successive  transformations  of  an  arbitrary  vector 
& by  means  of  the  matrix  D is  based  upon  the  fundamental  theorem* 

For  several  of  the  simplest  cases  this  method  has  been  investigated 
in  the  works  of  Kises  and  of  Duncan  cited  above,  and  — with  a more 
general  view  — in  the  x/ork  of  Aitken*  The  results  obtained  in 
56-10  of  the  present  work  are  partially  contained  in  Aitken's  work. 
^See,  e*g*,W,  p*  l£7* 
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6.  THE  GENERAL  CASE 
Let  us  consider  the  sequence  of  vectors 

a,  Da,  D2a,  ...  , D?*,  ...  (6.1) 

where  a is  an  arbitrary  vector.  If  r be  the  maximum  nunber  of 
linearly  independent  vectors  in  this  sequence,  the  vectors 
a,  Da,  •••  , Dar~^  will  be  independent,  and  having  deteri.iined  the 
coefficients  of  the  relation  of  linear  dependence  holding  ariong 
the  first  r + 1 vectors  of  the  sequence  (6.1) 

Y<prZ  + /■1Dr_1a  + •••  + jr'r.pa  + Yr-  * 0 <6«2) 

it  is  possible  to  write  the  characteristic  equation 

*0^  + ^l^1  + **•  + Yr-l*  + *r  ’ 0 (^3) 

for  the  matrix  defining  the  tr an s f or. i ati on  induced  by  transforma- 
tion L (with  matrix  D)  in  the  minimal  invariant  manifold  contain- 
ing the  vector  a.  T .is  is  the  substance  of  the  met  .od  elaborated 
by  A.  N®  Krylov  and  K.  N.  Luzin  for  composing  the  characteristic 
equation. 

From  the  preceding  paragraph  it  follows  that,  to  the  assigned 
degree  of  accuracy,  the  nunber  (which  we  shall  designate  by  m)  of 
successive  linearly  independent  vectors  in  sequence  (6.1),  begin- 
ning with  D^a,  will  diminish  with  the  growth  of  p.  If  the  coef- 
ficients of  the  linear  dependence  holding  among  the  vectors 
D^a,  9 EP+m§  be  determined: 

c DP+ma  ♦ c,DP'Hn”^a  ♦ •••  ♦ c ,D^+^a  ♦ c ■ 0 
o-l  - m-1  ~ m 3 


(6.U) 
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then  with  a sufficiently  large  p all  roots  of  the  equation 
m-1 


» m ji 

c0A.  + ClA 


+ c , A + C e=0 
m-1  m 


(6.5) 


will  be  equal  in  modulus  and  will  coincide  with  the  roots  of  equa- 
tion (6.3)  that  are  of  greatest  modulus. 

Talcing  p sufficiently  large,  the  degree  of  equation  (6.3)  may 
be  reduced  to  the  number  of  dominant  latent  roots  ’ aving  elementary 
divisors  of  highest  degree. 

Actually,  it  follows  from  formula  (U.U)  that  among  the  vectors 
D^i  corresponding  to  the  latent  roots  with  equal  moduli,  with  in- 
creasing p that  particular  vector  increases  most  rapidly  for  which 
the  degree  of  the  corresponding  elementary  divisor  is  greatest. 

Thus  having  taken  an  arbitrary  vector  a,  after  a sufficiently 
large  number  of  transformations  of  it  by  means  of  matrix  D,  by  the 
aforesaid  method  an  equation  (6.5)  can  be  obtained,  the  roots  of 
which  will  be  the  dominant  latent  roots. 

The  question,  how  numerous  must  the  transformations  be?  is 
difficult  to  decide  without  knowing  at  least  approximately  the 
magnitude  of  the  ratios  of  the  latent  roots  and  the  multiplicity 
of  the  corresponding  elementary  divisors.  Having  completed  a 
transformation  and  found  the  coefficients  of  the  relation  of 
linear  dependence  between  the  vectors  EPa,  •••  , Dp+ma,  the  ques- 
tion whether  p is  sufficiently  large  can  be  answered  by  testing 
whether  the  moduli  of  the  roots  of  equation  (6.3)  are  equal.  The 
necessary  and  sufficient  conditions  to  be  imposed  upon  the  coef- 
ficients of  the  equation  in  order  that  there  be  equality  of  moduli 
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of  the  roots  may  be  met  by  the  introduction  of  a new  unknown 

1 

jx  * | cQ  : c^j  A and  by  the  application  of  the  criterion  of  Cohn 

[9] 

for  the  distribution  of  the  roots  on  the  unit  circumference  • 
For  m « 2,  3,  U,  5,  a summary  of  these  conditions  is  set  forth 
below;  these  may  by  elementary  methods  be  obtained  readily  for 
the  given  most  simple  cases  (for  simplicity  cQ  has  been  taken 
equal  to  l).  In  contradistinction  to  the  general  criterion  of 
Cohn,  a preliminary  removal  of  the  small  roots  is  not  here  pre- 
supposed. 


1) 


2 

2.  Equation  A^  + c^A  + 


Condition  for  the  equality  of  the  moduli  of  the  roots: 


3 2 

2)  n = 3*  Equation  AJ  + c^<*  + C£  A + c^ 


0. 


Conditions  for  the  equality  of  the  moduli  of  the  rooi 


°3C1  » 


cx2  - c2l  < 2|c2I 


U 3 2 

3)  m • U.  Equation  A + c^A-5  + c^A  + A + • 0. 

Conditions  for  the  equality  of  the  moduli  of  the  roots 

c 2 

^ A v«i  A 

c-*c,  V^cT,  c^^lSc.,,  2 1 c..  1 - — 4 c 0 * +— — for  c.  > 0 

3141  3*  1 v^;  2 i*  u 


A ^ 

Cg  * 0,  c^  ■ — c^V-Cj^  — h V “C|^  * 0 for  * 0 


3 


1; 


.)  i.i  = >•  Equation 


+ °1 


Put  *»  ir  t 0.  Then  the  conditions  for  the  equality  of  the  no  luli 
of  the  roots  take  the  form 


c3=rc2,  r c19 


te-i 

2 

h-i 

r 

7 

r 

If  a check  reveals  the  condition  for  the  equality  of  t e moduli  of 
the  roots  for  equation  (6*5)  to  have  been  satisfied,  the  roots  nay 
then  easily  be  determined.  First  the  real  roots,  equal  to  -lcnlm» 
may  be  removed.  After  t’  is  separation  of  the  real,  the  complex 
roots  cr  then  be  determined  by  the  met  ods  ordinarily  applied  for 
the  solution  of  equations  by  the  method  of  Graeffe. 

Following  the  determination  of  the  latent  roots,  the  cor- 
responding invariant  manifolds  are  found.  Since  the  vectors  be- 
longing to  the  simplest  invariant  manifolds  satisfy  equations  of 
form  (U#6)  or  (U*12),  constituting  determinate  linear  combinations 
of  the  vectors  D^a,  •••  , D^+ma,  it  is  possible  to  obtain  the  vec- 
tors forming  the  basis  of  the  simplest  invariant  manifold  corres- 
ponding to  one  or  another  latent  root,  or,  more  correctly,  of  that 
part  of  the  manifold  that  is  contained  in  the  invariant  manifold 
defined  by  the  vectors  D^a,  •••  , D^+ma.  Actually,  letting  A* 
be  some  root  of  equation  (6.5)  of  multiplicity  k,  the  vectors 

(D  - A*)kDP+1a,  i - 0,  1,  ...  , m - k 

will  not  have  components  belonging  to  the  invariant  manifolds  cor- 
responding  to  A * Of  the  system  of  vectors  obtained,  the  com- 
ponents that  belong  to  the  invariant  manifolds  corresponding  to 
whichsoever  other  latent  roots  of  those  obtained  may  be  annihilated 
in  like  manner# 
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7.  so:  ;e  freu~:t  cases 


In  practice  cases  in  which  a significant  number  of  latent 
roots  coincide  in  modulus  are  rarolj'  encountered,  nor  are  cases 
in  which  elementary  divisors  of  high  degree  are  present.  A more 
detailed  investigation  of  the  process  of  determining  the  latent 
roots  of  the  invariant  manifolds  is  made  below  for  the  simplest 
cases.  This  will  be  sufficient  in  the  surpassing  majority  of 
practical  problems. 

1.  Simplest  will  be  the  case  where  A^  is  a simple  latent 
root  and  I / I In  this  case  for  sufficiently  large  p the 

direction  of  the  vector  D?a  will  coincide  with  the  invariant 
direction  corresponding  to  A^,  and  accordingly,  beginning  with  a 
certain  p,  the  direction  of  the  vectors  D^+1a  will  coincide  the 
one  with  its  predecessor  — i.e.,  their  coordinates  will  be  pro- 
portional. This  proportionality  of  coordinates  ser/es  also  as  an 
indicator  that  the  direction  of  the  vectors  of  our  sequence  EP+1a 
coincides,  within  the  limits  of  accuracy  of  the  computation,  with 
the  invariant  direction  of  the  matrix  D. 

The  latent  root  corresponding  to  t is  invariant  direction,  will 
be  equal  to  the  ratio  of  corresponding  coordinates  of  the  vectors 
DP+1a  and  l^a. 


For  computational  convenience,  in  place  of  the  sequence  of 
vectors  D^a  (p  * 1,  2,  •••),  a sequence  of  vectors  b^  could  be 
considered,  defined  by  the  recurrence  formula  Vi  ‘ P>DV  -o*a* 
Here  ^ are  arbitrary  numbers  selected  in  answer  to  the  d • . , for 
simplification  of  the  confutation.  With  the  growth  of  p,  b will 
approximate  in  direction;  the  ratio  of  correspond!-  coordinate! 
of  the  vectors  b^+^  and  b^  will  approach  ^ 


I 
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The  introduction  of  auxiliary  multipliers  should  be  recommended 
in  all  the  cases  following  as  well,  but  we  shall  not  do  this  out  of 
concern  for  simplicity  of  e::position. 

Vie  shall  obtain  just  such  a picture  if  A^  is  a multiple 
latent  root  of  multiplicity  r,  but  the  corresponding  elementary 
divisors  are  all  of  the  first  degree  (here  as  before  | \^\  fL  \ A^|  )• 
The  invariant  directions  obtained  in  the  manner  indicated  above 
will  here  ordinarily  be  distinct  for  distinct  initial  vectors  a; 
they  will,  however,  belong  to  one  and  the  sane  r-dinensional  mani- 
fold corresponding  to  A^,  in  which  any  direction  is  invariant 
(v.  § 2,  case  3). 

If  we  find  the  invariant  direction,  using  any  initial  vector 
whatsoever,  but  only  one  sue’",  then  those  two  cases  will  be  indis- 
tinguishable. 

Let  us  now  assume  that  with  a real  latent  root  At  there  cor- 

X X' 

respond  the  elementary  divisors  (A  - A^)  •••  , (A  - A^) 


*1**2 


* t; 


Then,  as  was  pointed  out  in  § 2,  the  component  of  an  arbitrary 
vector  a with  respect  to  an  invariant,  s^-dimensional 

m X^  + %2  + •••  + t^)  manifold  Eg^  corresponding  to  Ap  if 
only  it  does  not  equal  zero,  will  always  belong  to  a certain  axial 
invariant  manifold  lying  within  Eg^,  anc*  ^ * X^»  If  only 
| X 2I  / | A.^|,  then  on  the  strength  of  the  fundamental  theorem  it 
may  be  considered  that  the  vector  qQ  ■ E^a  for  sufficiently  large 
p,  and  within  the  limits  of  admissible  accuracy,  belongs  to  E^, 


. 
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The  equality  of  the  moduli  of  the  rcotc  of  equation  (6.5)  betokens, 
as  has  been  indicated  above,  the  fulfilment  of  the  condition  of  the 
condition  of  the  fundamental  theorem  and  the  fact  that  the  process 
of  computing  the  vectors  Dma  may  be  terminated.  In  the  case  in 
hand  equation  (6.5)  must  take  the  form 

( A - Xp1*  - 0 (7.1) 

which  may  be  easily  established  after  finding  the  maximum  number  k 
of  linearly  independent  vectors  in  the  sequence  qQ  * D^a, 
q^  » D^qQ,  (i  * 1,  2,  •••)  and  obtaining  the  equation  of  linear 
dependence  (6.U). 

Practically,  here  as  in  the  rest  of  the  simplest  cases,  it 
is  not  necessary  to  find  all  the  coefficients  of  linear  dependence 
(6.U).  Equations  (U.8)  and  (U.9)  are  consequences  of  equation 
(U.7),  and  the  two  former  may  serve  for  the  determination  of 
on  the  supposition  t’at  equation  (6.5)  reduces  to  (7«1)»  having 
determined  by  means  of  any  two  of  these  equations  (only  two 
or  three  of  the  coordinates  of  the  polyvectors  are  usually  sat*- 
isfied  in  the  determination),  one  can  judge  the  correctness  of 
the  assumption  made  by  the  coincidence  of  the  roots  obtained, 
however  this  coincidence  cannot,  strictly  speaking,  serve  as  a 
precise  criterion  whether  equation  (6.U)  has  the  form  (U.7),  and 
a definitive  verification  of  the  value  of  that  is  found  is  to 
be  made  simultaneously  with  the  verification  of  the  invariant 
direction,  discussed  below. 

Let  it  now  be  established  that  qQ  belongs  to  a certain  in- 
variant axial  manifold  E^  corresponding  to  the  latent  root  A^. 
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The  vectors  q^,  ^9  •••  , q^  may  serve  as  a basis  for  this 
manifold,  but  they  will  not  form  a fundamental  basis. 

Let  us  now  construct  the  groups  of  vectors  qi  according  to 
the  following  formulae : 


Si* 


(i  - 1,  2,  •••  , k) 


„ r „r—l  _ „r-l  / i * 1,  2,  •••  , k - r\ 

3i  * 3i+l  - *l9i  * \ r » 1,  2,  ...  , k - 1) 


(7.2) 


Since  the  vectors  satisfy  the  equation  (D  - jqjt}0  « 0, 

P 

the  vectors  q^  will  satisfy  the  equation 

<D - >i>k"V  - 0 


Moreover,  the  vectors  q.r  for  fixed  r will  be  linearly  inde- 
pendent  on  the  strength  of  the  linear  independence  of  the  vectors 
q^  and  may  accordingly  be  taken  as  a basis  of  the  axial  manifold 
In  particular,  the  direction  of  vector  will  coincide 

with  that  of  the  latent  vector  corresponding  to  A^*  The  vectors 

k-1  k— 2 1 o 

q1  , i i form  a basis  in  that  is  fundamental. 

In  fact,  they  are  linearly  independent  on  the  strength  of  the  in- 
dependence of  q^,  •••  , q^  and  satisfy  relations  (2*1),  which  fol- 
lows from  (7«2)*  The  computation  of  gives  conclusive 

verification  of  the  accuracy  of  the  computations  and  of  the  justice 
of  the  application  of  formulae  (4*8)  and  (4*9)  since,  obviously, 
only  in  the  case  of  the  correctness  of  all  assumptions  made  will 
we  obtain  = ^l9l"^* 


2k 


2 • Let  us  pass  now  to  the  next  case,  where  ^ * a 2 and 
l A^l  ^ I A^J.  Let  the  dominant  latent  roots  be  the  complex 
quantities  A-^  * and  A 2 - the  corresponding  elementary 

divisors  being  all  of  the  first  degree,  and  | \^\  / f • 

Then,  on  the  strength  of  the  fundamental  theorem,  with  suf- 
ficiently large  p the  vector  D^a  * qQ  will  belong  to  the  two- 
dimensional  imaginary-axial  invariant  plane  that  corresponds  to 
our  latent  roots,  and  equations  (U.13)  will  be  satisfied  accord- 
ingly. From  them  p and  0 may  be  determined.  For  basic  vectors 
there  may  be  taken 

f 1 * - f cos  0qo  + q1  and  m f sin  0% 


It  may  be  readily  verified  that  the  first  two  equations  of 
(3*1)  will  be  satisfied. 

If  the  complex  latent  roots  are  multiple  with  elementary 
divisors  (A  - A-^1,  •••  , (A  - 

(A  - A 2)^,  •••  , ( A - A2)^11,  the  vector  D^a  = q^  with  suf- 
ficiently large  p will  belong  to  a certain  imaginary-axial  in- 
variant manifold  where  X * Xm  depends  on  the  choice  of  a and, 
accordingly,  equations  (U.17 ) will  be  satisfied.  The  first  two, 
or  the  first  and  last  of  them,  give  the  possibility  of  readily 
determining  ^ and  0.  For  control  it  is  better  to  make  use  of  at 
least  three  of  the  equations. 

Vectors  ^2*  ***  * ^20  are  ^nearly  independent  and  can 
serve  as  basis  for  E ^ y Let  us  construct  now  the  vectors 
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Si0  - qA  (i  - 1,  2,  •••  , 2%) 

- (D2-2 f cos  jZSD  + <-2)qJ-1-qj7|-2r  cos  ♦ f2^1  (7.3) 

/r  - X,  2,  ...  , % - 1 \ 

V1  - 1,  2,  •••  , 2«  - r)J 

The  vectors  q^r[i  - 1,  2,  ...  , 2(X  ~ r)]  can  serve  as  a basis 
for  the  imaginary- axial  invariant  manifold  E2(^-r)  containe<i  in 

since  they  are  linearly  independent  and  satisfy  the  equation 

2 2 Y— r 

(D  - 2^»  cos  0D  + f a qi  » 0 which  follows  from  their  defini- 
tion, The  system  of  vectors  q^1*,  q2r,  (r  « %9  % - 1,  •••  , 2,  1), 
despite  that  it  will  not  be  a fundamental  basis  — since  equations 
(3,1)  will  not  be  satisfied  — will  nevertheless  possess,  as  may 
be  readily  confirmed,  the  most  important  property  of  a fundamental 
basis,  namely,  that  the  first  2k  vectors  of  this  system  farm  the 
basis  of  the  imaginary- axial  manifold  E2k  contained  in  E^, 

3,  Let  us  suppose  that  the  matrix  D has  several  distinct 
latent  roots,  each  of  these  having  the  same  modulus,  and  this  the 
greatest.  If  a limitation  to  real  latent  roots  only  be  made,  a 
case  is  here  possible  in  which  A^  ■ - A.2#  Let  the  corresponding 
elementary  divisors  be 

•••  , (A-A^)  m and  (A^A^)1*1*  •••  > (A+A^)^ 

tti  •••  4 Xm;  ^ <■  h2  4 ...  * h^ 
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We  introduce  the  designations:  ♦ •••  ♦ % ■ s.,  and 

± c m x 

hl  + h2  + *##  + \ " s2*  T}ien  ^he  components  a^  and  ag  the 
arbitrary  vector  a in  the  corresponding  invariant  manifolds 
and  Eg^  will  belong  to  certain  invariant  axial  manifolds  E^  and 
E^,  where  %2  4 For  definiteness  let  % ^ 0n 

the  strength  of  the  fundamental  theorem  the  vector  qQ  * D^a  with 
sufficiently  large  p will  satisfy  the  equation 


(D 


t,  u 

- Xx)  1(D  + X^)  ^ - 0 


(7.U) 


and  accordingly 


t,  t?t* 

(-1)  Ki  [5o9i  •'*  3^2-1 1 * [3i32  ***  (7,5) 


However  this  equality  as  well  as  others  analogous  that  may  be 
obtained  from  (7#U)  are  little  suitable  for  the  determination  of 
A. ^ and  the  fundamental  vectors. 

From  equation  (7*U)  it  follows  that  q^  and  q^  satisfy  the 
equation 


(D2  - X^)\  - 0 


(7.6) 


whence  it  issues  that 


[3oS2"*92i*.2S2r+2  **•  9efJ  “(^ 

(r  - 0,  1,  , %1-  1)  (7.7) 

©2  «,-*■) 

Ai  [3i33,*,32^-il 
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and  the  coordinates  of  the  polyvectors  [ •••  q^^gl  3111(1 
[qiq3  •••  n ] will  not  be  proportional.  Equations  (7*7 ) per- 

fectly coincide  in  essence  with  equalities  (U.8),  since  they  are 
obtained  from  (7»6)  just  as  (U.8)  is  from  (U.6).  From  them  A^ 
is  readily  found.  Having  determined  the  number  of  linearly  inde- 
pendent vectors  among  q^,  qg,  •••  # qg^  , which  number  must  be 
equated  to  X1  ♦ X2f  t2  also  can  be  determined.  For  the  determina- 
tion of  the  fundamental  bases  in  E^  and  let  us  proceed  in  the 
following  manner.  Construct  the  two  series  of  vectors 

f - £ - si  pi  - fi;i-  4~ht 


,) 


(7.8) 


Si  - Si7i+  xl\ 


Since 


if 


(D  - Ax)  x - 0 

h £ 

x lies  in  Ey  , and  p^  • (D  - A^)  q^f  which  follows  directly 
flrom  formulae  (7«7),  then  p.  lies  in  Ey  • In  the  same  fashion 
q^^  may  easily  be  shown  to  lie  in  ^ • 

The  further  determination  of  the  fundamental  bases  defining 
E^  and  E^  proceeds  Just  as  in  case  2 already  considered.  We  find 
the  vectors 

tft  tft- 1 % tj#*  , „ VN 

Pi  ■ pSi  ♦ *-iPi  “ x»  2»  •••  » 
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v v *+k -1 

With  % m %2  * 1 5111  the  vectors  win  coincide  in  direc- 

tion, which  will  also  be  the  invariant  direction  belonging  to  Ey  • 
For  the  fundamental  basis  the  vectors  (%  ■ 0,1,2,  “l) 

may  be  taken.  The  fundamental  basis  for  Ey  is  determined  anal- 

H2 

ogously  in  accordance  with  the  vectors 

It  is  to  be  remarked  that  had  the  number  been  determined 
in  the  manner  indicated  above,  it  would  have  been  sufficient  to 
take  only  ^ ^ ^he  vectors  q^°  and  p^°.  We  would  then  have 

obtained  the  vectors  p^*  defining  E^,  and  of  the  vectors 

defining 

In  conclusion  be  it  remarked  that  in  the  simplest  case,  when 
X^  • %2  m 1,  formulae  (7*7)  will  give  us 


52  " ^i25o*  53  " *i25i 


(7.9) 


The  invariant  directions  corresponding  to  and 
X^  ■ • X^  will  coincide  with  q2  ♦ X^_  311(1  ^ “ ^l9l* 

§8.  SEVERAL  FREQUENT  CASES  FOR  WHICH  THE  DEGREE 
OF  EQUATION  (6.5)  DOES  NOT  EXCEED  U 
Of  the  possible  simplest  cases,  we  shall  yet  consider  those 
for  which  the  degree  of  equation  (6#5)  does  not  exceed  four.  This 
is  possible,  excepting  those  cases  already  considered,  if 


1) 

xx  - X 2 • f 

X3  - 

I X <•  f,  where  t - ± 1 

2) 

X i ■ X 2 • 

^3  ■ f * 

3) 

X^«X.2"^e  » 

X3-  SU 

| x^l  < r 
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along  with  which  the  elementary  divisors  corresponding  to  the 
indicated  latent  roots  will  be  of  the  first  degree1. 

On  the  strength  of  the  fundamental  theorem,  with  sufficiently 
large  p the  vector  D^a  ■ will  satisfy  the  equation  (the  enumera^- 
tion  is  with  respect  to  the  cases  above) 

1)  (D2  - 2f  cos  0 + (»2)(D  - 0 , 

2)  (D2  - cos  0 ♦ f2)(D2  -f2)^  • 0 (8.1) 

3)  (D2  - 2f  cos  01+f2)(D2  - 2f  cos  02  + f 2)qo  • 0 


From  these  equations  the  fundamental  relations  permitting 
the  determination  of  /?  and  0 may  at  once  be  obtained: 


1)  tSiVfa1  ■ 4r3[3oSiS2]»  [3oS233]  " fr[3o3i335» 

[qoSi^l  * + 2 cos 


(8.2) 


2) 


“ 0 ^but  ^SoSiS^^0^  ^iq3Sl^ *^^t5oSiq3^ 

(8.3) 

[3o33a*J  “f2[3o3i9ul»  [303i3ii]  ’ 2f  cos  ^SoSii^ 


Case  3 may  be  considered  as  inciting  two  possibilities  omitted 
from  consideration,  to  wit:  the  presence  of  elementary  divisors 
of  the  second  degree  for  A3  in  case  1,  or  the  presence  of  elemen- 
tary divisors  of  the  second  degree  for  a single  pair  of  complex 
latent  roots.  Actually  the  corresponding  equation  for  ^ will 
coincide  with  equation  (8.1)  with  fa  * 0 or  fa  * TT  or  • #2* 

The  obtaining  of  these  values  of  fa  from  equations  (8.U)  will 
serve  as  a signal  of  the  presence  of  elementary  divisors  of  the 
second  degree,  in  accordance  with  which  corresponding  changes  in 
the  process  of  finding  the  invariant  manifolds  must  be  introduced* 
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3)  [2l32?39U]  ■rUt3o3il2,i3]»  t3oS2S33u]‘(,,2t3o3iS2Si11> 

* 2^(cos  0X+  cos  * (8#U) 

2 

[qoSi^Sl^]  - 2f  (1  ♦ 2 cos  0X  cos  WSoW&J 

A coincidence  of  values  for  obtained  from  the  first  two  of 
the  equations  of  each  group  can  serve  as  a control  of  the  correct- 
ness of  the  hypothesis  that  the  conditions  of  the  fundamental 
theorem  have  already  been  fulfilled  and  that  equation  (6.U)  has 
in  fact  the  form  of  one  of  the  equations  (8,1). 

The  determination  of  the  invariant  manifolds  is  here  possible 

by  the  general  method  indicated  at  the  end  of  f6.  Thus  for 

2 2 

example  , in  case  1 the  vector  (D  - 2p  cos  0 + pc) q^  will  give  the 
invariant  direction  corresponding  to  Ay  and  the  vectors 
(D  - ” ^“3^Sl  wi^‘  determine  the  imaginary-axial 

manifold  corresponding  to  g*  Analogous  formulae  may  be 
readily  obtained  for  the  rest  of  the  cases  as  well® 

§ 9.  PRACTICAL  APPLICATIONS  OF  THE  METHOD  OF  ITERATION 
In  computations  the  determination  of  the  dominant  latent 
roots  of  the  given  matrix  D and  of  the  corresponding  invariant 
manifolds  are  effected  in  the  following  manner. 

An  initial  vector  a is  selected;  this  ordinarily  has  only 
one  coordinate  different  from  zero,  or  has  all  coordinates  equal- 
ing 1.  The  vectors  D“a  (m  ■ 1,  2,  3,**#  ) »re  computed1.  The 
nature  of  the  change  of  the  vectors  of  this  sequence  points  to 

^Or 


rather,  the  vectors  b^  • see  f 7 
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one  or  other  case  of  the  presence  of  distinct  roots  of  greatest 
modulus*  If  the  vectors  Draa,  with  a large  m,  change  direction 
insignificantly,  this  points  to  the  presence  of  an  elementary 
divisor  of  degree  higher  than  the  first  for  latent  root  A or 
to  the  presence  of  roots  of  close  modulus*  If  the  methods  of 
§ 7 do  not  give  satisfactory  results  the  iterations  must  be  con- 
tinued for  the  separation  of  the  unequal  roots* 

The  presence  of  complex  roots  ( f 7,  case  2)  may  be  revealed 
by  strong  oscillations  of  the  directions  of  the  vectors  if1  a, 
marked  by  oscillations  of  the  signs  of  these  vectors1  coordinates* 

If  however  the  vectors  If  a.  If*2*,  D^a,  ...  slowly  change  direc- 
tion, as  do  vectors  •••  , then  case  3 of  J7  holds. 

If  with  sufficiently  large  m one  cannot  succeed  in  identifying 
one  of  the  three  possible  simplest  cases  considered  in  57,  one 
is  obliged  to  search  for  the  least  number  of  linearly  independent 
vectors  in  the  sequence  I^a,  Dra+^a,  •••  and  to  determine  the  coef- 
ficients of  the  relation  of  linear  dependence  (6.U)  between  them. 

For  this  compute  here  and  there  a few  of  the  determinants  of  second, 
third  and  higher  orders  that  are  the  coordinates  of  the  correspond- 
ing polyvectors,  until  a proportionality  of  the  coordinates  of  the 
different  polyvectors  of  the  first  order  is  obtained*  If  the  number 
of  linearly  independent  vectors  does  not  exceed  four,  the  formulae 
of  § 8 may  be  utilized;  in  the  contrary  case  the  general  method 
( f 6)^  must  be  employed. 

^See  the  examples  of  computations  at  the  end  of  this  work  ( f lU). 
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In  conclusion  we  make  one  more  remark  of  a practical  character. 
Instead  of  carrying  on  the  computation  with  a sequence  of  vectors 

a.  Da,  D2a,  ...  , D®a  (9.1) 

one  could  consider  the  sequence  of  matrices 

D,  D2,  DU,  D8,  ...  (9.2) 

by  means  of  which  analogous  results  would  be  obtained.  In  the 
simplest  case,  in  which  all  the  latent  roots  are  real  and  distinct, 
in  the  limit  there  would  be  obtained  a degenerating  matrix  the 
rows  of  which  would  be  proportional  to  the  coordinates  of  the 
latent  vector  corresponding  to  the  latent  root  At  first  sight 

the  second  method  would  seem  to  be  preferable  in  practice,  since 
it  secures  a more  rapid  convergence.  However  this  is  not  entirely 
true,  because  to  obtain  one  member  of  sequence  (9.2)  one  must  carry 
through  the  same  computational  work  as  one  would  to  obtain  n members 
of  sequence  (9*1) • 

Since  it  is  important  for  us  to  expend  the  least  labor  to  ar- 
rive at  the  power  m necessary  in  the  given  case,  obviously  the 
first  method  may  prove  to  be  the  more  advantageous. 

For  example,  if  the  degeneration  indicated  above  is  arrived 
at,  within  the  limits  of  computational  accuracy,  with  m - 16,  the 
first  method  would  be  more  advantageous  than  the  second  for  n > U; 
if  it  would  be  sufficient  to  take  m • 30,  the  first  method  would 
be  the  more  advantageous  for  n > 6, 

A combination  of  the  computational  methods  is,  however,  most 

2 U 

expedient.  At  first  find  D , and  perhaps  even  D , and  then  consider 


’ 
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the  sequence 

a,  D2a,  D^a,  •••  or  a,  D^a,  D8a,  Dl6a,  •••  (9.3) 

2 1 

The  invariant  manifolds  of  the  matrices  D and  D coincide  , 

2 

and  the  latent  roots  of  D are  equal  to  the  square  of  the  latent 
roots  of  matrix  D.  Therefore  to  obtain  results  in  sequence  (9*3) 
one  can  take  half  as  many  members  as  in  (9.1),  which  appreciably 
reduces  the  quantity  of  computations. 

To  assist  visualization,  a table  is  subjoined  exhibiting  the 
quantity  of  computations  for  various  n and  m by  diverse  methods. 

For  a unit  there  has  been  adopted  the  quantity  of  computation  neces- 
sary to  complete  one  transformation  of  the  vector  by  means  of  the 
matrix  Dj  (obviously  these  units  are  not  of  equal  value  for  dif- 
fering n). 


m • l6 

m - 32 

m ■ 61* 

Computational 

Method 

n*l* 

n * 6 

n ■ 8 

n -1* 

n *6 

n * 8 

n*  1* 

n*  6 

n ■ 8 

Iteration  by  D 

16 

16 

16 

32 

32 

32 

61* 

61* 

61* 

■ « D2 

12 

11* 

16 

20 

22 

21* 

1*0 

1*2 

1*1* 

it  « 

12 

16 

20 

16 

20 

21* 

21* 

28 

32 

« « D8 

1 1* 

20 

26 

16 

22 

28 

20 

26 

32 

Succeeding 

Involution 

16 

21* 

32 

20 

30 

1*0 

21* 

36 

1*8 

if  D has  latent  roots  + X and  -A.  and  the  corresponding  invariant 
directions  f^  and  f2,  then  for  D2  all  vectors  *1^  ♦ will  be 
latent  vectors  corresponding  to  the  latent  root  In  an  anal- 

ogous fashion  the  invariant  manifolds  will  be  united  in  the  case 
for  which  + X and  - X are  multiple  latent  roots.  However,  when  the 
invariant  manifolds  of  D and  the  corresponding  latent  roots  are 
known,  the  determination  of  them  for  D represents  no  particular 
difficulty.  (See  | ll*.  example  1). 


f 10  PROPERTIES  OF  THE  INVARIANT  MANIFOLDS 
OF  CONJUGATE  MATRICES 

By  using  the  methods  indicated  in  the  preceding  paragraphs 
one  can  determine  the  dominant  latent  roots  and  the  invariant 
manifolds  constructed  on  axial  and  imaginary-axial  manifolds  cor- 
responding to  the  latent  roots  found.  At  the  same  time,  depending 
on  the  selection  of  the  initial  vector  a,  there  may  remain  un- 
discovered individual  latent  roots  with  the  same  modulus,  or  the 
multiplicity  of  the  latent  roots  found  will  not  be  fully  determined. 
The  latter  will  always  occur  if  the  elementary  divisors  of  the 
matrix  D - A E are  not  relatively  prime. 

Let  the  invariant  manifold  found  have  s dimensions;  the 
aggregate  of  basic  vectors  for  the  simplest  axial  and  imaginary- 
axial  manifolds  entering  Eg  forms  a basis  for  Eg.  We  shall  desig- 
nate them  by  fp  fg,  •••  , fg.  We  shall  show  how  to  pass  from  a 
matrix  D to  a new  n - s-dimensiomal  matrix  D^,  defining  a linear 
transformation  in  E g and  possessing  the  same  latent  roots  as  Df 
with  the  exception  of  those  already  found  (taking  into  account  the 
multiplicity). 

To  accomplish  this  we  shall  prove  an  auxiliary  theorem  con- 
cerning the  matrix  Dc,  the  conjugate  of  D [i,e.,  transpose  - 
translator].  The  elementary  divisors  of  the  matrix  Dc  - A E will 
obviously  coincide  with  the  elementary  divisors  of  D - AE,  and 
accordingly  the  structure  of  the  invariant  manifolds  of  the 

matrices  D and  D will  be  identical.  The  connection  between  these 
c 

manifolds  is  given  by  the  following  theorem. 
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Any  axial  invariant  manifold  of  the  matrix  Dc  corresponding 
to  the  latent  root  is  orthogonal  to  all  the  invariant  manifolds 
of  the  matrix  D,  with  the  exception  of  those  which  contain  the 
axial  manifolds  corresponding  to  that  particular  latent  root.  Any 
nr-diraensional  axial  invariant  manifold  of  the  matrix  Dc  is  orthog- 
onal to  any  k - m-dimensional  invariant  axial  manifold  of  the 
matrix  D that  lies  within  whatsoever  k-dimensional  axial  manifold 
corresponding  to  that  particular  latent  root. 

We  shall  first  prove  the  first  half  of  the  theorem.  For  this 
it  is  sufficient  to  prove  that  the  vectors  g^,  g2,  •••  , g^  that 
form  the  fundamental  basis  of  an  arbitrary  axial  manifold  of  the 
matrix  D corresponding  to  the  latent  root  X . are  orthogonal  to 
the  basic  vectors  h^,  h2,  •••  , of  an  arbitrary  axial  manifold 
of  the  matrix  D corresponding  to  a latent  root  i The 

basis  hp  h2,  •••  , we  shall  consider  as  being  fundamental. 

From 

h * DcSi  “ xAl*Si  ' fi  * Dh  * Vgi  * -i 

it  follows  that  g^  • h^  • 0.  Furthermore,  from 
Jl * + #il  * Dcll  " *i^##i  W-2,3,**«,k) 

there  follow  sequentially 

gx  • - Ji  * J 3 " * Jl*i?k  " 0 

In  exactly  the  same  way  we  find 

h' h mh.'hm  *•*  ■ 2?i*  J»  * 0 
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Whence  by  the  method  of  complete  induction,  on  the  basis  of  the 
equality 

V Dcfr-Ir  * D£*  ’ V ' Jh-1 '^-1  * Jr  * 0 

it  follows  that 

gy*  h^  ■ 0 (r*  1,2, •••,m,  l,2,»*»,k)  (10*1) 

which  is  what  was  required  to  be  proved* 

The  second  part  of  the  theorem  is  proved  analogously*  Let 

gl*  g2*  9 ^ 311(1  Si9  2 9 ***  9 now  desiSnate  the  actors 

forming  the  fundamental  basis  of  the  axial  invariant  manifolds 
of  the  matrices  Dc  and  D corresponding  to  the  same  latent  root  X • 
From 


II  * D£i  * fl  * ^ /i-1^  " * Ddl  " *Ji  * Jl 

it  at  once  follows  that 

g^*  ^ • 0 (i  • 1,  2,  •••  , k - 1) 

Furthermore,  on  the  ground  that 

§2  * D£i  ■ Jh  * ( XIi+ii-l)  “ Ji  * Dof2  ■ h * ( *#2  + Jl} 

we  obtain 


g2#/i  ” 0 (i  - 1,  2,  •••  , k - 2) 
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Proceeding  further,  we  will  find  that  for  any  X * m 

gj«  f±  - 0 (i  . 1,  2,  •••  , k - X)  (10.2) 

which  proves  the  second  half  of  the  theorem. 

The  proved  theorem  holds  both  for  real  and  complex  space,  as 
is  seen  immediately  from  its  proof. 

Let  now  f^,  fg,  , f^  and  g^,  j2,  •••  , g^  be  the  fundamental 
system  of  generators  for  a certain  maximum  axial  manifold  of  the 
matrix  D and  for  the  corresponding  manifold  E’^  of  the  matrix  Dc* 

By  maximum  we  here  understand  such  an  axial  invariant  manifold  as 
cannot  be  contained  in  an  axial  manifold  of  a greater  number  of 
dimensions. 

The  product  g^«  f^  will  be  different  from  zero.  Let  us  assume 
at  first  that  the  elementary  divisor  ( X - to  which  the  invar- 

iant manifolds  correspond  is  relatively  prime  to  all  the  rest  of 
the  elementary  divisors.  Then  if  it  be  assumed  that  g^  • f^  * 0, 
the  vector  would  be  orthogonal  to  n independent  directions, 
which  could  not  be,  as  g^  / 0, 

If  now  for  the  latent  root  we  have  r elementary  divisors, 
given  the  vectors  in  hand  f^,  •••  > the  vectors 

§1,  gg*  •••  , g^  will  not  be  uniquely  determined.  In  particular, 
for  the  vector  g^  there  may  be  taken  any  vector  of  some  r-dlmen- 
sional  manifold,  and  this  selection  made  so  that  f 0, 

which  follows  from  considerations  analogous  to  those  presented 


above. 
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It  is  also  easily  shown  that  • f^^  ^ 0 (i  ■ 1,2, 1). 
Assuming  the  contrary,  in  the  simplest  case  we  arrive  immediately 
at  the  conclusion  that  gi+^  may  be  linearly  expressed  in  terms  of 
gp  •••  , g^,  which  is  impossible  on  the  ground  of  the  linear  in- 
dependence of  the  vectors  g^,  •••  , g^.  In  a case  where  the  vector 
g^^  is  not  uniquely  determined,  it  may  always  be  chosen  such  that 
the  demanded  inequality  is  satisfied* 

From  the  preceeding  it  follows  that  it  is  always  possible  to 
choose  numbers  c^  such  that  the  vectors 

i 

li  - t °ikik  (i  • x»  2»  (10-3) 

k*l 

will  form  a system  of  vectors  reciprocal  to  f^,  •••  , f^,  that  is, 
such  that 

Ii+1*  (i,k  • 0,1,  ...  , *-l)  (10.U) 

Finding  such  a system  of  reciprocal  vectors  in  each  simplest 
invariant  manifold  of  the  matrix  Dc,  we  shall  obtain  a complete 
system  of  vectors  reciprocal  to  the  system  of  n vectors  composed 
from  the  fundamental  basis  of  each  simplest  invariant  manifold  of 
the  matrix  D. 

In  the  case  of  complex  latent  roots,  the  fulfilment  of  equa- 
tion (10.U)  for  the  complex  vectors  determining  the  corresponding 
invariant  manifolds  leads  to  the  fulfilment  of  analogous  equations 
for  the  real  vectors  forming  the  fundamental  basis  of  the  cor- 
responding imaginary-axial  manifolds. 
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S 11.  THE  TRANSFORMATION  OF  THE  GIVEN  MATRIX  D 

Let  us  designate  by  g^',  , gg*  the  vectors  corresponding 

to  fg,  •••  , fg,  selected  from  the  system  of  reciprocal  vectors 
constructed  in  the  way  indicated  above. 

Furthermore,  we  designate  by  (j  * 1,  2,  •••  , n)  the 

coordinates  of  the  vector  f^,  and  by  p ^ (j  * 1,  2,  •••  , n)  the 
coordinates  of  the  vector  j^'» 

We  introduce  into  the  space  new  coordinate  vectors,  namely^" 


~V  * ~6i  - P\h  - <>2% 

~ n-s+k  * 


“ /^s^s  (i*l,2,**-,n-s) 

(11.1) 

(k«  l,2,3,»»*,s) 


This  substitution  of  coordinates  is  equivalent  to  the  trans- 
formation of  the  matrix  D by  means  of  a matrix  T by  the  formula 

D*  - T*1  DT  (11.2) 

Here  T denotes  a matrix  whose  elements  t^  are  defined  by 
the  following  equations: 

hk  " Sik  " ^ <*%  t1"1*2.  •••  » "S  k-1,2,  ...  , n-s) 

(11.3) 

h.,  n-s+k  ’ (i-1.2,  •••  , n;  k- 1,2,  ...  , s) 

We  designate  the  elements  of  the  matrix  by 

^Here  it  is  assumed,  without  detriment  to  the  generality,  that  the 
determinant  of  the  s-th  order,  | ©<  J"S4’*C| , is  different  from  zero. 


ho 


On  the  basis  of  formulae  (11.3)  we  readily  find 


rUc-sik  (i.k-1,2,  •••  ,»-.) 


^i,n-s+k  " * 


,i  - V 

(i-1,2,  ••• 

» nj  k-1. 

2,  • ••  , 

.n-s+l 

. n-s+k-1  i 

. n-s+k+1 

_ n 

• • • 

*1  *1 

• • • 

^ n-s+l 

_,n-s+k— 1 . i 

n-s+k+1 

_ n 

*<2  *2 

*2 

*'*«2 

.n-s+l 

/m  A a 

^n-s+k-1  ^ i 

n-s+k+1 

...  11 

s 

s s 

*s 

s 

. n-s+l 
*i 

^ n-s+2 

^ ••• 

M n 
*1 

. n-s+l 
**2 

n-s+2 

•x.  2 ••• 

n 

* 2 

n-s+l 

n-s+2 

n 

*s 

*8 

(11.U) 


The  elements  of  the  matrix  D*,  transformed  on  the  basis  of 
(11*2) , may  be  computed  by  the  formula 

n n 


E Z t 


d±k'  'A  ,ti 


(11.5) 


Substituting  in  the  expressions  exhibited  above  for 

and  r v,  we  arrive  at  the  following  results: 

(*K 


15  V 


u 


n 

Z T,  d v (i,j  »l,2,»».,n- s)  (11.6) 
*- n-s+l 


2)  d^i  - 0 (i«l,2,«»»fn- s;  J«n-s+l,***fn)  (11.7) 

3)  The  s-dimensional  matrix  with  elements  d^*(i,J  > n-s) 
will  be  the  matrix,  in  Jordan  canonical  form,  of  the  transformation 
induced  in  the  invariant  manifold  Efl. 
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U)  If  some  axial  manifold  Ey  present  in  Eg  and  corresponding 
to  the  latent  root  A.  with  basis  f -,»••••  f is  not  maximum  — 
that  is,  there  exists  outside  of  Eg  a vector  f such  that 
Df  • Xf  + — then  the  element  d^ * for  i ■ n-s+m+r  and 

j n-s  will  coincide  with  the  j-th  contravariant  coordinate 
(with  respect  to  the  new  coordinate  vectors)  of  the  vector  g cor- 
responding to  f in  the  reciprocal  system  of  vectors* 

For  other  values  i * n - s and  j < n — s,  d^.*  * 0. 

Points  2,  3,  and  h may  even  be  demonstrated  without  computa- 
tion by  formula  (11,5)#  and  directly  on  the  strength  of  the 
properties  of  the  fundamental  bases  of  the  simplest  invariant 
manifolds  and  of  formula  (11*1)  with  the  coordinate  vectors 
replaced. 


§12.  THE  INVARIANT  MANIFOLDS  OF  THE  TRANSFORMED  MATRIX, 

Let  us  consider  the  n-s  ^dimensional  manifold  E constructed 

n— s 

on  the  vectors  •••  • X ' • This  manifold  does  not  have  common 

1 7 7 'n-s 

directions  with  E ; accordingly  together  with  E it  completely  d^* 
s s 

fines  the  whole  space.  The  matrix  with  elements  d^*  (i,jin-s) 

defines  a transformation  induced  by  matrix  D in  E . Any  vector 

n-s 

belonging  to  E has  its  last  s coordinates  (in  E ) equal  to  zero, 
n-s  n 

Its  first  n-s  coordinates  are  transformed  by  the  matrix  just 

as  they  are  by  matrix  D‘ , Consequently  any  axial  manifold  E^  of 

the  matrix  D not  contained  in  E„  will  lie  in  E _ and  will  there 

s n-s 

be  an  invariant  axial  manifold  of  the  matrix  D^, 

If  a certain  axial  manifold  Er  belonging  to  Eg  is  part  of  an 
axial  manifold  of  a greater  number  of  dimensions,  E^#  the  matrix 

will  have  the  invariant  axial  manifold  corresponding  to 

the  same  latent  root. 


Indeed,  let  h^,  •••  , be  a fundamental  basis  in  E^.  Then 
h^,  ...  9 hr,  under  the  condition,  will  coincide  with  certain  of 
the  vectors  f.<>  Let  h.  ■ f . (i  * 1,  2,  •••  . r)«  Consider  the 
vector  a * + c2hr+2  + • ••  + where  c p ...  , c^r  are 

arbitrary  numbers#  Then 

“ XS  * °2-r+l  + c3~r+2  4 ***  4 %r^-l  4 CA  (12*1) 

In  the  new  system  of  coordinates  the  last  s coordinates  of  the 
vector  a are  equal  to  zero,  since  h^^,  ...  , h^s  according  to  the 
condition,  lie  in  . The  vector  h has  all  coordinates  equal  to 
zero  except  the  (n-  s ♦ m ♦ r)-th.  Since  the  matrix  D transforms 
the  first  n*»  s coordinates  just  as  does  D*  , it  thence  follows  that 

DA  • 4 cAm  4 — 4 %A/-1  (12-2) 

where  a#,  h#r+^,  • • • $ h*^_i  signify  the  n - s -dimensional  vectors 
the  coordinates  of  which  coincide  with  the  first  n-  s coordinates 
of  the  vectors  a,  ...  , h^^.  It  follows  from  equation  (12.2) 

that  the  vectors  h^r+^,  •••  , are  a fundamental  basis  of  the 
axial  manifold  E^_r  of  the  matrix  D#.  The  vector  h^r+1  is  the  axis 
of  this  manifold. 

f 13.  FURTHER  DETERMINATION  OF  THE  LATENT  ROOTS 
AND  THE  INVARIANT  MANIFOLDS  OF  THE  MATRIX  D 
If  a certain  invariant  manifold  Eg  and  the  corresponding  latent 
roots  have  been  determined,  one  can  go  over  to  a matrix  of  n - s 
dimensions  to  reveal  the  remaining  latent  roots  and  invariant  mani- 


fold E , of  the  matrix  D . This  may  be  done  using  the  very  method 
S'  c 


' 


U3 


by  which  Eg  was  found.  We  only  remark  that  to  make  possible  the 
future  determination  of  the  vectors  g^  it  is  necessary  that  the 
axial  invariant  manifolds  entering  into  Eg|  be  certainly  maximum. 

It  is  always  possible  to  attain  this  by  taking  as  the  initial 
vector  during  the  determination  of  Eg|  a linear  combination  of 
vectors  directed  along  the  invariant  directions  already  found,  ly- 
ing in  E 5 (given  the  presence  in  E of  imaginary-axial  manifolds- 
the  first  two  generating  vectors  of  the  fundamental  system  are  to 
be  taken  from  them).  Then  on  the  strength  of  the  results  obtained 
in  f ID,  this  initial  vector  will  have  components  different  from 
zero,  lying  in  the  maximum  axial  manifolds  of  the  matrix  D^,  which 
manifolds  correspond  to  those  axial  manifolds  of  the  matrix  D from 
which  we  have  taken  the  axial  directions.  Therefore,  generally 
speaking,  s'  * s. 

Further  transformations  are  to  be  conducted  as  indicated  in 

S10,  11. 

During  this,  since  for  Dc  the  dominant  latent  roots  are  known, 
as  soon  as  the  vectors  Dcma,  iT1*,  •••  , are  computed,  the  vectors 
forming  the  fundamental  basis  for  Eg|  may  be  found  at  once  by  the 
corresponding  formulae. 

We  now  remark  that  in  formula  (11*6),  by  which  the  elements 
of  the  matrix  D#  are  computed,  the  coordinates  of  the  vectors  g^ 
do  not  appear.  Its  [D#]  latent  roots  and  invariant  manifolds  may 
accordingly  be  found  even  without  the  computation  of  the  vectors  g^. 

However  to  obtain  the  invariant  manifolds  of  the  matrix  T one 
must  return  to  the  original  coordinate  system,  which  cannot  be  done 
unless  one  has  the  coordinates  of  the  vectors  g^. 


■ 


In  this  manner,  after  the  determination,  by  the  method  indicated 
in  the  first  half  of  the  present  work,  of  the  invariant  manifold  Eg 
of  the  matrix  D,  which  manifold  is  composed  from  the  axial  or  imagi- 
nary-axial simplest  manifolds  corresponding  to  the  dominant  latent 
roots,  we  construct  the  matrix  with  elements  which  are 

defined  by  formula  (11*6) • The  latent  roots  of  this  n- s -dimen- 
sional matrix  coincide  with  n-  s (with  due  regard  for  multiplicity) 
as  yet  unfound  latent  roots  of  the  matrix  D*  By  the  same  method 
the  dominant  latent  roots  of  the  matrix  and  the  corresponding 
invariant  manifolds  can  be  determined. 

If  we  want  to  obtain  these  manifolds  in  the  original  system 
of  coordinates,  we  must  carry  through  the  process  indicated  in 
f 12  for  determining  the  vectors  g^» . To  obtain  the  coordinates 
of  any  vector  in  the  original  system  of  coordinates  the  formula 

n-6 

\ • Z tit  O'  - 1.  2>  — , n)  (13.1) 

will  obviously  serve.  Here  a#^  signifies  the  coordinates  of  the 
corresponding  vector  in  and  t^  are  determined  by  formula 

(11*3)  • 

Continuing  in  this  manner,  we  shall  pass  to  matrices  of  ever 
lower  order,  gradually  determining  all  the  latent  roots  of  the 
matrix  D* 

iJi.  NUMERICAL  EXAMPLES 

Example  1.  Let  it  be  required  to  determine  the  latent  roots 
and  invariant  directions  of  the  matrix 


U5 


U — 6 0 — U 

-6  U h 0 

0 U U 6 

-U  0 6 U 


For  expediting  the  convergence  we  shall  carry  out  the  iters^ 
2 

tions  not  with  D but  with  D : 


68 

-U8 

-U8 

-32 

-U8 

68 

32 

U8 

-U8 

32 

68 

U8 

-32 

U8 

U8 

68 

We  shall  take  for  the  initial  vector  a « {l,  0,  0,  oj  • 

In  Table  1 are  given  the  vectors  * “a;  the  computations 

terminate  at  k ■ 8,  since  the  ratio  of  the  coordinates  of  the 
vectors  ag  and  coincide.  The  magnitude  of  this  ratio  gives 

Xj2  • 196.000. 


Table  1 


k 

he 

he 

1 

1 

68 

-U8 

-1*8 

-32 

2 

0.1 

978 

-960 

-960 

-896 

3 

0.01 

1873 

-1860 

-1860 

-181*1* 

k 

0.1 

361*93 

— 361*1*2 

-361*1*2 

-36389 

5 

0.01 

71W*1* 

-71U25 

-711*25 

-711*07 

6 

0*01 

11*0000 

-139993 

-139993 

-139987 

7 

0.01 

271*389 

-27U387 

-271*387 

-271*381* 

8 

0.01 

537799 

-£37798 

-537798 

-537797 

A direct  check  immediately  convinces  that  f^  (1,  -1,  —1,  -1) 
defines  the  invariant  direction  of  D corresponding  to  the  latent 
root  * Hu 


For  the  determination  of  the  remaining  latent  roots  we  pass 


to  three-dimensional  space;  in  this  case  formulae  (11*6)  will 
yield 

d.Uc-dik-^d!*  (M-l,  2,  3) 


where  tXy  •<  ^ are  the  coordinates  of  f.  Performing  the 

computation,  we  obtain 


D - 


* 


36 

-16 

-16 


0 

0 

20 

-16 

-16 

20 

Taking  for  the  initial  vector  b • jl,  0,  o]  , we  compute 
b^  • b (Table  2)  up  to  k • 6,  when  the  ratio  of  coordinates 
gives  * 36.000. 


Table  2 


k 

he 

4 

1 

1 

36 

*16 

-16 

2 

1 

1296 

-6U0 

— 6I4O 

3 

0,1 

U666 

-2330 

-2330 

h 

0,1 

16798 

-8398 

-8398 

5 

0.1 

60U73 

-30236 

-30236 

6 

0.1 

217703 

-108851 

-108851 

For  the  coordinates  of  the  latent  vector  there  can  be  adopted 


1 

'*2 


O • f 


2 

‘*2 


•<-  O ■ ■ J I 


l*2 


J 


The  presence  of  two  zeros  in  the  first  row  of  at  once 

o 

indicates  that  for  D the  latent  vector  corresponding  to  • 36 

*■  C 6 


, 


will  be  g#c  (1,  0,  0).  Passing  on  to  a matrix  of  the  second  order 
by  formulae  (11.6),  we  obtain 


2 

Direct  computations  will  give  the  latent  roots  • 36  and 

2 

• U,  and  the  corresponding  invariant  vectors  ( 1,  -1)  and 

W1* 0)* 

In  order  to  obtain  the  invariant  directions  of  the  given  matrix 
one  must  perform  a reverse  transition  to  the  space  of  the  greatest 
number  of  dimensions.  First  we  find  the  invariant  directions  D# 
using  formula  (13*1): 

X i 

* *k*=  ±J^  2 ^i^-wtk  e ^ " 1>2#3) 

here  t#^_  » ^ **2  2*  fi*2  are  the  coordinates  the 

vector  normalized  so  that  £+2*  2 * 

The  computations  give 


*11  * *12  “ *32  * 0 * *21  * *31  ’ °*S  * *22  “ 1 


whence  we  find 

f3  - {o, -0.5.  0.5}  5 fu-  (0,0.5,  0.5} 


For  the  further  transition  to  four-dimensional  space  we  are 
served  by  the  formulae! 

- z **  <k  - t « i,2,3,u) 

K i-1,2,3  M 


. 


U8 


where  ®<  and  is  the  coordinate  of  f^, 

which  must  be  normalized  so  that  f^  • f^  ■ 1,  Performing  the 
computations,  for  convenience  replacing  f^by  2f^  ^ and  f^  by  Uf^, 
we  obtain: 


■11" 

0.75 

"12 

K 

0.25 

"l3 

- 0.25 

•21  " 

0.25 

"22 

■ 

0.75 

*23 

« -0.25 

•31" 

0.2S 

"32 

« 

—0.25 

"33 

- 0.75 

'Ul  * 

0.25 

V 

° 

-0.25 

"U3 

- -0.25 

0, 

-li  » 

~f3" 

r -1,  1,  0] 

l. 

Ik  * I1. 1. 

For  control  and  the  determination  of  the  signs  of  A^#  we  com- 
pute Df^: 

Df2  - {o,  6,  -6,  o}  , Df3  - [6,  0,  0,  6 } , Df ^ ^ 2,  2,  2,  -2}  . 

From  the  results  it  is  obvious  that  • 2 and  that  the  invar- 
2 

iant  plane  for  D , each  direction  of  which  plane  is  invariant,  con- 
tains only  two  invariant  directions  corresponding  to  the  latent 
roots  A 2 * 6 and  A^  - - 6;  (see  the  last  footnote  to  f?)»  As 
the  vectors  defining  these  invariant  directions  we  may  take 

Df2  + 6f2  ■ (—6,  6,  -6,  -6 1 , Df2  “ 6f2  ■ ^6,  6,  —6,  6 j • 


Example  2»  Let  it  be  required  to  determine  the  latent  roots 
and  invariant  manifolds  of  the  matrix 

r 


6 -3  U 1 

U 2 1*  0 

h -2  3 1 

U 2 3 1 


h9 


We  shall  iterate  the  matrix 


uu 

-30 

27 

11 

U8 

-16 

36 

8 

32 

-20 

20 

8 

U8 

-12 

36 

8 

D2  - - 


having  adopted  for  the  initial  vector  a jl,  0,  0,  oj-  • 
Sk  * hA  are  set  forth  in  Table  3. 


The  vectors 


Table  3 


k 

Pk 

1 

1 

Uh 

U8 

32 

1*8 

2 

1 

1888 

2880 

1U72 

3072 

3 

1 

70208 

122112 

56832 

133632 

k 

0.01 

2U302 

U5312 

120101 

50196 

5 

1 

80U811 

1566708 

675012 

171*7956 

6 

1 

2586328U 

518U7680 

21903680 

58111*512 

7 

0.1 

81321309 

166532333 

69366118 

1872711*05 

The  vectors  a^  turn  out  to  be  nearly  of  the  same  direction, 
but  they  do  not  tend  so  rapidly  to  the  limiting  invariant  direction 
as  in  the  preceeding  example,  the  which  compels  us  to  surmise  that 
multiple  roots  are  present  (case  1 of  those  considered  in  $ ?)• 
Actually  the  coordinates  of  the  polyvectors  [a^a^],  [a^a^],  [s^a^] 
prove  to  be  approximately  proportional  and  their  ratio  gives,  by 
formulae  (U«8)  and  (Lu9)  for  % m 2 


— - - 13.7081  , 2 A 2 - - 


2 


5U.8323 


50 


p 

whence  A * 27*Ul62,  and  the  values  obtained  from  both  formulae 
coincide  sufficiently  well  to  indicate  the  correctness  of  the  suj>- 
mise  made  about  the  multiplicity  of  the  root#  [See  Note  3 at  the 
end  of  the  translation.] 

2 

The  values  of  Aj  obtained  may  be  easily  verified  and  made 

more  precise.  Actually  the  vector  q « (D^  - X^)a£  ought  to  give 

the  invariant  direction.  The  computation  shows  that  the  ratio  of 
2 

the  coordinates  of  D q and  q is  in  fact  constant,  but  equal  to 

27.U166.  If  the  arithmetic  mean  of  both  values  obtained  — 27.U16U  ♦ 

2 2 

be  taken  for  A,  , and  q and  D q again  be  found,  we  will  obtain  for 

the  ratio  of  coordinates  27.U16U2.  Again  taking  the  arithmetic 

2 

mean  27#Ul6Ul  — - for  A-^  , we  shall  obtain  for  the  ratio  the 
quantity  27 «1±16U06 • Limiting  ourselves  to  five  figures  of  decimals, 

*»  27 *416141  may  be  adopted.  For  the  vector  f^(*^,  > *1^ 

defining  the  invariant  direction,  there  nay  be  taken  the  vector  q,  the 
first  coordinate  of  which  has  been  reduced  to  unity;  and  for  the 
vector  £2^*  2^ 9 *2  9 *2  9 *2^  9 c*e^iniri2  together  with  the  axial 
invariant  manifold  E2>  the  vector  a^,  normalized  likewise: 

f1  - (1.0000000,  2.3U16U11,  0.89UU272,  2.6832819) 

f2  • (1.0000000,  2.0U78311a,  0.8529882,  2.3028576) 

Formulae  (11.6)  for  the  transition  to  two-dimensional  space 
reduce,  in  the  case  in  hand,  to  the  form 


Ti*d.y 


(i,J  - V) 
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Here 


(-D 


k-1 


tlk" 


(k  • 3,U;  i • 1,2) 


The  computations  give 


. |“°- 

l-1* 


1*60365 

01*3965 


1.01*3957 

1.627556 


The  characteristic  equation  for  will  be 


X - 1.167191  A + 0.31*0585 


both  roots  of  which  are  equal  to  0.583596.  The  invariant  direction 
of  the  matrix  for  this  latent  root  is  defined  by  the  vector 
(1.01*3961,  1.01*3957 ) , which  can  with  sufficient  accuracy  be  re- 

placed by  the  vector  f#^(l,  1).  Any  other  vector  linearly  inde- 
pendent of  even  f^j^O,  1)*  defines  together  with  the 

corresponding  axial  manifold. 

To  obtain  the  vectors  and  f^  in  foujvdimensional  space. 


the  invariant  manifolds  of  the  matrix  D must  be  found, 

c 


The 


computations  b^  • are  given  in  Table  1*.  Far  the  initial 

vector  b the  vector  (1,  2,  1,  2),  approximately  coincident  in 
direction  with  has  been  taken.  (See  above.) 
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Table  h 


k 

Pk 

1 

1 

268 

-106 

191 

51 

2 

1 

1526U 

-10776 

9076 

U036 

3 

1 

638528 

-515U56 

391008 

186592 

k 

1 

235U2016 

-20167808 

121*21312 

7200960 

5 

0.01 

8109220 

- 718U133 

1+1725U1 

25U5979 

6 

0.1 

26769560 

-2U23330U 

135U2622 

85U7652 

7 

0.1 

85831325 

-78677820 

1*2891*709 

27732092 

8 

0.1 

269U21boU 

-250357391 

1331*09375 

878135U2 

o 

The  invariant  direction  of  is  defined  by  the  vector 
- (T?  - A^)by,  and  for  the  vector  defining  the  corresponding 
axial  invariant  manifold,  b,  » h^  may  be  taken.  From  these  vectors 
we  construct  the  vectors 


Jh//*1  9 ^1  * 9 ^1^  and  ^2^2  9 P 2 9 9 P 2 ^ * 


which  are  reciprocal  with  f^  and  f^,  by  the  formulae 


h 5x(52*I2)  *2 

|2  - * #i  - - Ttp^TTIvT?  + 


The  computations  give 

g1  - (-6.2611961*,  6.7611970,  -2.6766166,  -2,3013866)  , 

g,  - ( 6.981*7961*,  -6.987961*,  3.237630 9,  2.1*131876)  . 

Further,  we  find  t^  by  formulae  (11.3)  and  the  coordinates 
of  f^  and  fj^  by  formulae  (13.1): 
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tn  - 0.276U00,  t21  - 0.357789,  t31  - -0.35776U, 
t12  - 0.223600,  t22— 0.528611,  t32  - -0.089U5, 
f3  = (0.500000,  -0.170822,  -O.UU721U, 

f^  - (0.223600,  -0.528611,  -0.089U50, 


- 0.71556U 

- -2,057206 
-1.3U16U2) 
-2.057206) 


A control  computation  D£3  gives  for  the  ratio  of  the  coordinates 
of  Df3  and  f the  numbers:  0.76393,  0,76396,  0.76393,  0.7639U, 
which  are  in  full  agreement  with  the  value  found  earlier  for 


x32  - 0.583596. 

Example  3«  To  be  determined,  the  invariant  directions  and 
latent  roots  of  the  matrix 


D • < 


15 

1 

7 

7 

17 


11 

3 

6 

7 

12 


- 9 

- 3 

- 3 

- 3 

-10 


-15 

- 8 
-11 
-11 
-16 


The  computations  a^  «=  for  si  c (1*  1*  1,  1,  1)  are  given 
in  Table  5»  The  multiplier  is  • 0.1  for  k ■ lU  — 16,  and  • 1 
for  all  the  other  values  of  k. 


. 


Table  5 


k 

2k 

0 

1 

1 

1 

1 

1 

1 

8 

2 

5 

5 

8 

2 

7 

-20 

-5 

-8 

7 

3 

-178 

-130 

-154 

-169 

-178 

4 

-833 

-23 

-485 

-461 

-818 

5 

761 

2660 

1502 

1964 

836 

6 

191*71 

9679 

15211 

16369 

19352 

7 

52211* 

-18508 

23669 

18137 

49904 

8 

-190157 

-21*3932 

-206891 

-249068 

-195947 

9 

-1596136 

-1*69192 

-1133416 

-1170457 

-1568476 

ID 

-I842395 

2651*723 

- 23993 

- 640231 

-1631510 

11 

223326U1 

17637221* 

20113532 

-22992248 

22517846 

12 

101*98231*9 

711*6589 

66159973 

63693665 

101661229 

13 

- 1*7763868 

-2821*78951* 

-134608699 

-193622083 

- 62157448 

1 4 

-19561*1621* 

-102855319 

-157227482 

-172014507 

-194403670 

15 

-51*521511*6 

152017762 

-275505683 

-221133520 

-515708456 

16 

156676258 

222035520 

177876001 

22062834 6 

164069771 

The  course  of  the  computations  points  to  the  presence  of  complex 
roots,  and  the  lack  of  proportionality  of  the  coordinates  of  the  bi- 
vectors points  to  an  equality  of  the  moduli  of  more  than 

two  latent  roots*  For  determining  them  we  shall  employ  the  foimulae 
of  5 8.  Having  taken  qQ  ■ a^,,  we  compute  one  after  another  the 
coordinates  of  the  polyvectors  * ##*  > 

t301iS2l»  tSiS2S3}»  •"  • 

Then  as  the  ratios  of  the  coordinates  of  the  bivectors  coincide, 
within  the  limits  of  the  first  two  or  three  significant  figures,  the 
trivectors  turn  out  to  be  noi>-proportional.  Finally,  two  coordinates 
of  four  different  vectors  are  computed;  these  are  determinants  of  the 


■ 


fourth  order,  composed  of  the  first,  second,  third,  fourth,  and  of 
the  second,  third,  fourth  and  fifth  coordinates  of  the  correspond- 
ing vectors.  The  results  of  the  computations  are  brought  together 
in  Table  6* 


Table  6 


Polyvector 

First 

Coordinate 

Second 

Coordinate 

1 90*32%! 

1.8881*21  xio22 

-2.360593  XlO22 

[3o3iS23!*1 

1.133201 

-1.U16U86 

[9o%.%9u! 

7.365237 

-9.20661*0 

[SoS2S39u1 

1.699807 

-2.12U73U 

[9i9?939u1 

U.2U8808 

-5.311186 

The  ratios  of  the  coordinates  coincide  sufficiently  well; 

(it  must  be  borne  in  mind  that  a loss  of  accuracy  occurs  during 
the  computation  of  the  determinants,  owing  to  the  subtraction  of 
numbers  little  different). 

On  the  basis  of  formulae  (8.U)  we  find  (the  means  of  the  two 
values  are  everywhere  adopted) : 


’[3oSiS2S3]’ 


22U.99U5 


[3o3A9u] 


15.000 


2f(co6«L+coS^.}^|j-6f)0068J  39.0018 


whence 


f-3.87297,  cos  J^-cos  02-O.3873Uf  - 1.50016  ♦ 13  ^706U 
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The  vectors  forming  the  fundamental  basis  of  the  correspond- 
ing  imaginary- axial  manifold  may  readily  be  found,  as  was  shown 
at  the  end  of  f 9» 

The  fifth  latent  root  may  at  once  be  obtained  from  the  trace 
(spur)  of  the  matrix: 

s*=A1+X2+X3  + ^|i+X^*l5+3+6-6-l6e55x^»  - 1.00032 

The  corresponding  invariant  direction  will  be  defined  by  the 
vector^ 

2 

« d^-a^cos  j*D  + f2)  a - {13.30,  22.05,  19.09,  16.09,  28.12  } 

By  a direct  check  it  is  readily  confirmed  that  the  vector 
{l3,  22,  19,  16,  28}  exactly  defines  the  invariant  direction  cor- 
responding to  the  latent  root  X^  ■ - 1.00000* 

(Jieceived  at  the  editorial  office  October  16,  19U0) 


■^n  the  work  of  Aitken  [U],  whence  the  matrix  for  the  present 
example  has  been  taken,  the  vector  is  given  incorrectly,  with- 
out any  indication  whatsoever  of  the  manner  in  which  it  was  ob- 
tained. There  are  also  given  the  exact  value  ^ and 

cos  <t  -\^755  . 


$ 

' 


. 
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Notes 

by  editor  of  translation 

1.  For  a definition  and  discussion  of  polyvectors  (or  multi- 
vectors)  in  the  restricted  sense  used  here,  see  Elie  Cartan, 

La  Geonetrie  des  e spaces  de  Riemann,  Mem.  des  sci.  math.f  do.  9q 
Paris,  1925,  page  U.  We  may  summarize  the  definition  briefly  as 
follows: 

Given  the  k column  vectors  q^,  q2,  •••  , q^  in  n-dimensional 
Euclidean  space,  form  the  matrix  Q of  n rows  and  k columns,  whose 
i-th  column  consists  of  the  n components  of  q^  (i  ■ 1,  2,  •••  , k). 
The  polyvector  [q^,  q2,  •••  > qk]  is  a (£)-dimensional  array  whose 
components  are  the  k-dimensional  minor  determinants  of  Q.  The 
components  of  a multivector  are  not  necessarily  independent  func- 
tions of  the  vectors  q^,  , q^. 

One  may  easily  prove  the  following  properties  of  polyvectors, 
which  are  used  in  proving  (U.9); 

i)  If  any  q±  - 0,  [q1,  • ••  , qfc]  - 0. 

ii)  If  ary  two  q’ s are  proportional,  [q^,  •••  , q^]  * 0. 

iii)  [q^,  •••  , 2^3  is  a linear  function  of  each  single  one  of 
its  arguments  q^,  •••  , q^. 

When  k is  2 or  3,  the  author  uses  the  terns  bivector  and 
trivector,  respectively, 

2,  Probably  the  author  meant  to  say  that  the  ratio  will 

grow  when  | * | X^| , except  when  a^  coincides  with  the  invariant 


direction, 


■ 


■ 
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3.  Here  and  later  the  author  divides  one  polyvector  by 
another,  according  to  the  following  rule:  Given  two  proportional 
vectors  a,b  , with  b / 0 and  a - •<  b ( ■ scalar),  the  quotient 

a/b  is  defined  to  be  • 


O 
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